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Abstract: In this paper, we study interacting extended Chaplygin gas as dark matter 
and quintessence scalar field as dark energy with an effective A-term in Lyra manifold. As 
we know Chaplygin gas behaves as dark matter at the early universe while cosmological 
constant at the late time. Modified field equations are given and motivation of the phe¬ 
nomenological models discussed in details. Four different models based on the interaction 
term are investigated in this work. Then, we consider other models where Extended Chap¬ 
lygin gas and quintessence field play role of dark matter and dark energy respectively with 
two different forms of interaction between the extended Chaplygin gas and quintessence 
scalar field for both constant and varying A. Concerning to the mathematical hardness of 
the problems we discuss results numerically and graphically. Obtained results give us hope 
that proposed models can work as good models for the early universe with later stage of 
evolution containing accelerated expansion. 
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1 Introduction 

A huge effort of fundamental physics to explain the dynamics of the universe opens big 
window for different speculations. However, the most intriguing effect in the recent uni¬ 
verse can not be explained properly and gives different philosophical fight between different 
minds. If we accept that Einstein was correct with his general relativity theory to explain 
accelerated expansion of the universe could be explained by negative pressure working 
against gravity. The belief of Einstein to the static universe made him to think about 
negative pressure which will stop the attraction of the gravity. However, we know that 
we have non static universe, moreover we know that we have accelerated expansion. Ac¬ 
cording to the last observational data analysis we estimated the amount of the negative 
pressure in our universe, which we call it as dark energy. The simple question about the 
nature of the dark energy still one of the intriguing questions and left free space for new 
speculations. A big class of the dark energy models are introduced in life by hand and 
used with the field equations of general relativity. This is working but is not a satisfac¬ 
tory approach, therefore modification of the field equations on the Lagrangian level could 
be considered more fundamental and satisfactory. This approach opened a wide range of 
the modified general relativity, but still with a big uncertainty, because the way and the 
form of the modifications plays a crucial role in our understanding and depends on the 
form of the modification of the same class of the modified general relativity gives different 
forms of the dark energy. This is model dependent approach and still left a big room for 
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speculations, of course some of the approaches could be combined with particle physics 
and can give a hope that this link will give us possibility to illuminate phenomenology and 
have a fundamental theory. Consideration of the modified general relativity promises new 
insight into our understanding of the universe and one of the hot and modern subject for 
the study. In this work we will consider one of the modifications and would like to see 
the behavior of the universe. Our analysis is based on exotic type of fluid which is called 
Chaplygin gas (CG) and sued to describe accelerating expansion of universe [1, 2]. This 
is indeed based on CG equation of state and developed to the generalized Ghaplygin gas 
(GCG) [3]. It is also possible to study viscosity in GGG [4-9]. Then, GGG was extended 
to the modified Chaplygin gas (MCG) [10]. Recently, viscous MCG is also suggested and 
studied [11, 12]. A further extension of CG model is called modified cosmic Chaplygin gas 
(MCCG) which was proposed recently [13-15]. The MCG equation of state (EoS) has two 
parts, the first term gives an ordinary fluid obeying a linear barotropic EoS, and the second 
term relates pressure to some power of the the inverse of energy density. However, it is 
possible to consider barotropic fluid with quadratic EoS or even with higher orders EoS 
[16-18]. Therefore, it is interesting to extend MCG EoS which recovers at least barotropic 
fluid with quadratic EoS. This is called extended Chaplygin gas (ECG) [19-24]. 

Another interesting models to describe dark energy are quintessence scalar field [25-30], a 
phantom field [31-37], the combination of both fields in a unified scenario which is called 
quintom [38-40], or more complicated models such as in K-essence [41]. 

Here, we would to combine ECG and quintessence field as the following scenario. We as¬ 
sume EGG as dark matter which is reasonable assumption at least in the early universe, 
however with possibility to extension in present epoch. The other component, dark energy, 
described by quintessence field. An important point is possibility of interaction between 
dark matter and dark energy. Gosmological models including Chaplygin gas usually used 
for unification of dark matter and dark energy, therefore we have a section of single ECG 
which can play role of both dark matter and dark energy. 

On the other hand a unified treatment of cosmological models is given within the context of 
Lyra manifold [42]. In that case there are some recent works of cosmology including Lyra 
manifold [43, 44] where A considered as a variable. So, we construct a new cosmological 
model and confirm it using the recent observational data. Some advantages of this work 
are as follows. 

1- First of all we use recently proposed type of Chaplygin gas. 

2- Reducing n -|- 2 free parameters of the ECG model to 3. 

3- Interacting quintessence and EGG is a new model. 

4- Gomparison with the recent BIGEP2 data. 

5- Several types of interaction term. 

6- Presence of Lyra manifold. 

7- Varying A. 

The paper is organized as follows. In next section we write down field equations which 
govern our model. Then, in section 3 we introduce the ECG and represent some cosmolog¬ 
ical behavior of the single EGG model. In section 4 we introduce our main models based 
on different forms of interaction term between quintessence field and ECG. In section 5 we 
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have numerical analysis of the case of constant A which extend to the case of varying A in 
section 6. Finally in section 7 we summarize our results and give conclusion. 


2 Field equations 


Field equations that govern our model are given by [43], 

1 3 3 

Rfiu ~ ‘fa — (^-l) 

where is a characteristic of Lyra geometry. Considering the content of the universe to 
be a perfect fluid, we have. 


Ti,u = ip + P)Uf,Uy - P9^u, ( 2 . 2 ) 

where = (1,0,0,0) is a 4-velocity of the co-moving observer, satisfying Uf^u^ = 1. Let 
be a time-like vector field of displacement, 

99^= (^^/3,0,0,o) , (2.3) 

where /3 = /3(t) is a function of time alone, and the factor is substituted in order to 
simplify the writing of all the following equations. By using FRW metric for a flat universe, 

= —dt^ + a{t)‘^ + r^idO"^ + sin^ 6d(j)^)) , (2.4) 

field equations can be reduced to the following Friedmann equations, 

37/2 -/32 = p + A, (2.5) 

and, 

2i7 + 377^ + ;32 ^ _P + A, (2.6) 

where H = d/a is the Hubble parameter, and an over dot stands for differentiation with 
respect to cosmic time t, and a(t) represents the scale factor. 

The continuity equation reads as, 

p + A + 2/3/3 + 3H(p + P + 2/3^) = 0. (2.7) 


With an assumption that, 

p + 3H{p + P) = 0. 

The equation (2.7) will give a link between A and f3 of the following form, 

A + 2/3/3 + 6 P /32 = 0. 


( 2 . 8 ) 

(2.9) 


In order to introduce an interaction between the dark energy and dark matter, we should 
mathematically split the equation (2.8) into the following equations, 

Pdm P^HipDM + Pdm) = Q, (2-10) 


and. 


PDF + 3H{p£)E + Pdf) = -Q, 


( 2 . 11 ) 


where Q is interaction term and expressed explicitly in the section 4. For the Dark energy 
sector we assume quintessence field which also introduce in the section 4. Before that we 
introduce the extended Chaplygin gas in the next section. 
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3 Extended Chaplygin gas 


In this section we review evolution of CG EoS which is given by the following primary 
relation, 

PcG = -, (3.1) 

PCG 

where B is an arbitrary parameter which usually considered as a constant. 

In order for the universe passes from a dust dominated epoch to a de Sitter phase through 
an intermediate phase described as mixture of cosmological constant and radiation the 
Chaplygin gas equation of state generalized to [45], 

Pggg = —, (3.2) 

Pgcg 

where a and B are free parameters. It is clear that a = 1 reproduces the pure Chaplygin 
gas model given by (3.1). This model is called the generalized Chaplygin gas model. There 
is also a class of equations of state that interpolate between standard fluids at high energy 
densities and generalized Chaplygin gas fluids at low energy densities which is called the 
modified Chaplygin gas with the following equation of state [46], 

Pmgg = ^Pmcg -X-, (3.3) 

Pmgg 

where A, a, and B are free parameters of the model. The case of ^ = 0 recovers generalized 
Chaplygin gas equation of state, and A = 0 together a = 1 recovers the original Chaplygin 
gas equation of state. The MCG equation of state has two parts, the first term gives an 
ordinary fluid obeying a linear barotropic equation of state, while there are some models 
with the quadratic equation of state [16-18], 

P = Po + u;ip + u)2P^, (3.4) 


where Pq, cui and W 2 are constants. Easily we can set Pq = a ;2 = 0 to recover linear 
barotorpic equation of state. 

Modified Chaplygin gas include only barotropic fluid with linear equation of state, while 
it is possible to extend them to including quadratic barotropic equation of state given by 
(3.4). It yields to introducing extended Chaplygin gas with the following equation of state, 

Pegg = ^ APecg -s— > (3-5) 

Pegg 

where a and B are constants. It is obvious that the n = 1 reduced to the modified 
Chaplygin gas with A = 1/2. Moreover, 


OJECG = 


Pegg 

Pegg 


^iPECG l+a ’ 
i=i Pegg 


(3.6) 


^2 _ dPEGG 
dpEGG ’ 


is EoS of ECG, and, 


(3.7) 









is sound speed. It is one of the important point of this paper which relate the coefficient 
of p%(jQ to n, so we reduced n free parameters of the model to one. 

In the next section we consider ECG as dark matter and assume quintessence field as dark 
energy to construct an interacting model. 


4 The models 


We assume quintessence scalar field as dark energy described by a scalar field (p and self 
interacting potential V{4>). Energy density and pressure given by, 


Pp 




(4.1) 


and, 

= (4.2) 

We consider models for the universe where an effective energy density and pressure given 

by, 

P = P<j> + PECG, (4.3) 


and. 


P = P^ + Pecg, 


(4.4) 


where p^j, and P^j, are energy density and pressure of the quintessence dark energy in the 
universe, while pegg and Pegg represent energy density and pressure of the ECG. An 
important point is possibility of interaction between components. 

The solving strategy and structure of the problem allow us to assume that the form of the 
potential V{<p) is given as follow. 


V ((/>) = Voexp 



(4.5) 


where A and 7 are arbitrary constants. Moreover, to make the system of the differential 
equation closed we assume two different situations, namely once we consider A to be a 
constant, while for the second case A is a dynamical quantity and we consider one of the 
classical form, 

A = p = P(j, + Pegg- (4.6) 

In this work we will consider three different types of interaction term [47] which gives us 
six different models depend on varying or constant A. At the first step we will perform 
numerical analysis of the important cosmological parameters in the case of constant A then 
consider varying A. Different forms of interesting interaction terms between quintessence 
and EGG provide us a huge number of the models. We start with the interaction term Q 
given by, 

Q = 3Hbp^, (4.7) 

then we consider the second model where Q reads as [48], 


Q = hHp^cg, 


(4.8) 
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where b, X and m are arbitrary constants. So, the parameter b can be interpreted as 
strength of interaction. The case of 6 = 0 gives non-interacting model. 

We start analysis of the models from the next section and start with the case where A is 
a constant 

5 Cases of constant A 

Consideration of the constant A can simplify field equations, and before graphical analysis 
we would like to provide mathematical forms of the simplified field equations. According 
to this assumption, the equation given by (2.7) will be modified as, 

p + 2p$ + 3H{p + P + 2/3‘^) = 0. (5.1) 

So, using the equation (2.8) we have, 

$ + 3HI3 = 0. (5.2) 

An integration of the last equation reveals dependence between f3 and the scale factor a as 
follow, 

15 = % (5.3) 

where uq is an integration constant. 

5.1 Model 1 

The first model has the interaction term given by the equation (4.7) with the constant A 
which yields to the following equations, 

1 , , Kl , ^HbV{^) ^ dV{^) ^ 

(j) + 3H{l + -)4>5 -^-^-= 

2 (j) d(p 

and, 

Pecg + 377(1 + ujecg)pecg - 3776(-(^^ + V{<f))) = 0, (5.5) 

where ujegg is the EoS parameter of the ECG given by (3.6). 

We solve the equation (5.4) numerically and fit the following function of the scalar field, 

</) = 2.1(1-(x-1)^), (5.6) 

Numerical result and best fitted value drawn in the Fig. 1. We can see that the scalar field 
is increasing function of time and yields to a constant value at the late time. 

Then, we solve the equation (5.5) numerically and find typical behavior of the extended 
Chaplygin gas density which illustrated in the Fig. 2. We can see interesting behavior, the 
energy density increases first to reach maximum value and then decrease to the infinitesimal 
value corresponding to the current value. 

We find good result corresponding to the case of n = 3 which is good news for our 
model which is a confirmation of the ECG, at least to the third order. Moreover, we can see 
far results away than observations for the case of 6 = 0 and n > 1, which means presence 
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Figure 1. Typical behavior of the scalar field against t (10^^ s) for the model 1. Dashed line 
represent behavior of the function given by the equation (5.6), while dotted line represent numerical 
solution of the equation (5.4). 





Figure 2. Typical behavior of the pecg against t (10^^ s) for the model 1. 
of interaction is necessary in the ECG model. 

Consideration of interaction terms bring us to investigate two important parameters = 
and ^ecg = Pecg/^H'^ to see solving coincident problem of the cosmological 
constant model. In our model, which cosmological constant is presence, we can check and 
find values of and ^egg at the same order. 

We also able to investigate evolution of the scalar potential of the model which presented 
by the Fig. 3. We can see that the scalar potential is decreasing function of time and 
vanishes at the late time as expected. 
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Figure 3. Typical behavior of the scalar potential against t (10^^ s) for the model 1. 

5.2 Model 2 

The second model obtained using the interaction term given by the equation (4.8) and 
gives two following differential equations, 

+ + = 0, (5.7) 

(p d(p 

and, 

PECO + 377(1 + loecg)pecg - = 0, (5.8) 

9 

which describe the dynamics of the content of the universe, uegg is the EoS parameter for 
the extended Chaplygin gas. Here, we examined higher order terms of the ECG and find 
good behavior of the Hubble and deceleration parameters. Analysis of total EoS parameter 
tells that ojtot ~1 at the late time for higher order terms of ECG. It means that Higher 
order terms are necessary. 

Moreover, we find that values of 91ecg and 91^ are of the same order. Finally we can obtain 
time evolution of the scalar field and scalar potential similar to the previous models. 

6 Cases of varying A 

Among the numerous possibilities concerning to the phenomenological form of A in this 
section we will consider one of the classical form given by the equation (4.6) proposed 
and considered for a long time. The analysis of the first model of this section with A 
starts by the interaction Q = SHbpq between two components. Mathematical formulation 
describing the dynamics of the energy densities of the components is the same as for the 
first two models, with one difference that the form of /3, therefore the form of the Hubble 
parameter, is different. In next subsections we will provide and discuss behavior of the 
cosmological parameters. 
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6.1 Model 3 


This model given by the equations (4.7) and (4.6). We find that the Hubble expansion 
parameter is decreasing function of time and yields to a constant at the late time which 
is in agreement with observations. Study of the deceleration parameter shoe acceleration 
to deceleration phase transition at the near past. The deceleration parameter begin with 
positive value (acceleration phase) and yields to minus one at the late time. Current value 
of this parameter q ~ —0.8 is available in this model. 

Also behavior of the EoS parameter is completely expected which yields to minus one at the 
late time. As before, we can check and find values of and ^ecg at the same order. Like 
to the previous models we find scalar potential and scalar field decreasing and increasing 
function of time respectively. 

Unfortunately, stability analysis of this model based on squared sound speed show that the 
model is completely unstable and it bring us to consider the next model. 

6.2 Model 4 

This model given by the equations (4.8) and (4.6). Behavior of cosmological parameters of 
this model is approximately similar to the previous model but we find that this model is 
stable at least before the late time. So, it can describe early and current state of universe. 
Having and give us scalar field and potential using the equations (4.1) and (4.2). 
So, hrst we should solve the equation (2.11) to find pecg then use it in the equation 
(2.10) to obtain scalar energy density. Therefore we can obtain evolution of /3, Hubble and 
deceleration parameters. 

The master equation is, 

Pegg + 377(1 + ojegg)pegg + H^p^cG = 0- (6.1) 

In the Fig. 4 we draw typical behavior of the ECG energy density for various values of m. 
We can see that the condition m > 1 is necessary. Also we find that evolution of pegg 
for higher m is faster than that with lower m. It is obvious that pegg behaves like some 
power of exponential so it help us to find typical behavior of (see Fig. 5). We plot p^ 
for m = 2 and find that higher m decreases initial values of density dramatically. However, 
at present we can say that and ^egg are at the same order. 

Now, using the relation, 

Cs = j> 0 , ( 6 . 2 ) 

where Cg is sound speed, it is easy to check stability of the model, we find that the model 
is stable at least from beginning to the late time. 

7 Conclusion 

In this paper we proposed a new cosmological model to describe evolution of universe. Our 
model is based on new proposed extended Chaplygin gas. In order to have a comprehensive 
model we considered the case of varying A in Lyra manifold. We studied also cases of 
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Figure 4. Typical behavior of the pecg against t (10^^ s) for the model 4. m = 0 (dash dot red), 
TO = 2 (solid black), to = 3 (dash green), to = 4 (dot blue). 



Figure 5. Typical behavior of the against t (10^^ s) with to = 2 for the model 4. 

constant A and compare our results with the case of varying A to find varying model 
as better case. We considered two special kind of interaction terms and concluded that 
the second kind yields to a good model. We found that an interacting ECG model with 
Q = H^bp^fjQ and A = p is an stable model which can describe evolution of the universe. 
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